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ABSTRACT 



In most behavioral science research very little attention is ever ' 
given to the probability of committing a Type II error, i.. > the 
probability of failing to reject a false null hypothesis. Recent 
publications by Cohen (1970, 1969) have led to a great deal of insight 
on this topic for the fixed-effects analysis of variance and covariance. 
It is the purpose of this presentation to provide social scientists 
with some insight in dealing with Type II error, and, therefore, 
optimum sample size and number of levels, in the random-effects 
analysis of variance. 



Introduction 

TLa random- effects analysis of variance Is described In most 
statistics texts (Winer, 1972; Glass & Stanley, 1970; Dixon & Massey, 
1969; Kirk, 1968; Hays, 1965; Guenther, 1964) used by social scientists. 
Each of these texts provides the reader with a procedure for selecting 
sample size for the fixed-effects model. However, none of them provide 
the reader with guide lines for selecting the number of levels (treatments) 
or elements within a level (subjects) for random-effects designs. 
Instead they provide examples la which it appears as though sample 
size and number of levels were chosen because they "looked good". 
It is exactly this kind of "looks good" procedure which these same 
authors have tried to avoid by discussing sample size selection for 
the f ixed-eff ebts model. 

Although the random-effects model is not used as often by social 
scientists as is the f ixed-ef f ectr model, it certainly does have a 
wide variety of applications. In this model a researcher is faced 
with the problem of having to draw Inferences about an entire set of 
distinct treatments or factor levels. In this case the researcher is 
not interested in the values of the individual treatment effects, 
as in fixed-effects designs, but in the variance of the population 
from which these effects were randomly selected. Possible populations 
which a social scientist might consider for study would be schools, 
teachers, psychologists ^ sociologists, time periods (e.g.. Glass & 
Stanley, 1970, pp» 452-462), animal strains, or classes. 

For example, suppose that a population of school teachers is 
available to teach reaUlng to first graders using a certain method. 



It is decided tha^: the method will adopted for use provided its success 
is not heavily dependent on the personalities of individual teachers. 
A random samplo of teachers is draxm and randomly assigned to a random 
sample of first grade pupils. The dependent variable (a standardized 
reading tes':) is selected and a one-way analysis of variance, random- 
effects mode},, is to be used to test the hypothesis th,Et there are 
no significant differences among teachers in the teacher population. 

Researchers faced with the preceding problem are often ignorant 
of how many teachers to select and the number of students to assign 
CO each teacher. The following discussion provides an answer to the^e 
questions for the one-way random-effedts analysis of variance, and 
indicates the problems involved in answering these questions for 
more donlplex random-effects designs. 

Background 

Cohen (1970, 1969) discusses the relationship between power (the 
probability of rejecting the null hypothesis (H^) when it is indeed 
false), sample size, and effect size (a measure of the effects one 
desires to detect) for the fixed-effects model. The calculation of 
power in that model is difficult because when Hq is false the F statistic 
is no longer distributed as a central F but as a noncentral F. In 
the random-effects model this is not the case. Here power calculations 
can be made using the central F distribution. 



Hypotheses 



la random effects designs, the null hypothesis is generally written 



with the alternative being 



A 0 



Here 9^^0. is a preassigned constant, and e is Cohen's (1969) "effect 
size" for random effects designs. That is, 6 is an index of the 
degree of departure from the null hypothesis which we want to detect. 



Parameter Selection 

The major problem faced by investigators is that of choosing 
values for 6 and 6^. These values may be based on the theory in the 
area which is being investigated. However, most often, they are informed 
hunches based on work done in a pilot study or on research found in 
the literature. 

Most beginning statistics texts discuss the rather limited case 
where 6^ = ^ This leads to the null hypothesis 

■Hj,: e.= 0. 

This hypothesis may be unrealistic since some differences most likely 
do exist among levels in the random-effects design. Therefore, values 
of of 0.10, 0.50 and 1.00 along with 0.00 were considered in this 
paper. 

It is the writer's firm belief that, if a researcher is unable to 
select values of e and 9^ then he should consider his experiment nothing 
more than a pilot study. Not knowing what 6 and 6^ probably are 
prevents one from being able to selact an appropriate number of observations 
with which to test Hp. Using too small of a sample, the researcher 
may decide to fail to reject a false null hypothesis. Type II error. 
Using too large of a sample, he may decide to reject a true null hypothesis. 
Type I error. In either case, he does not know enough about the area 
in which he is working to select an appropriate sample size. Unfortunately, 
if such errors occur they may find their way into print and lend 
confusion to an area of research. 



The One-Way Design 

In the one-way random-effects design, the null hypothesis is 
generally written 



with the alternative being 



2 



Here is the variance of the effects for factor A; is the variance 
of the sample elements, and 9o > 0 is a preassigned constant. 

Then the power of the F test (Scheffe. 1959) is a function of 
e = OA/ag. That is. the ratio ae^ae" is an index of the degree of 
departure from the null hypothesis which „e want to detect. So that. 
Power = Pr {F(L-1,N-L) > F(a; L-1, N-L)(l-h»e„)/(14ne) 

where 

F = the F statistic, 

L the number of levels, 

a = the level of significance, 

K » the total number of elements, 

n = N/L, the number of elements in a level of the design. 
The above formula for power is adapted from formula (7.2.12) in Scheffe, 
1959, p. 227. 

That is, poxver for any one-way random-effects analysis of 
variance is found using the central F distribution and is calculated 
by finding the probability of drawing an F value from a central 
F distribution, with L-1 and N-L degrees of freedom, that is greater 
than or equal to the value F(a;L-l,N-L) (l+neQ)/(i+ne) . 



6 

Therefore, given 9 and Qq, one may select the optimum number 
of lev6ls (L) to be used with a fixed total number of people (N) such 
that the power of the statistical F test is maximized. Here, N/L 
* n is the number of elements in each level of the design. Equal 
n's in each level have been assumed since some difficulties arise 
in the random-effects analysis unless there are equal numbers of 
observations (Hays, 1965, p. 419). 

Tables 1, 2, and 3 were constructed^ using the preceding 
procedure with the Central F distribution. In a computer program 
described in Barcikowski (1972, see Appendljl A). In these tables the 
optimum number of levels for the one-way random-effects model were 
found for given values of N, a, 9, and G^,. The procedure used was to 
select o, 8, 9^ and the total number of observations and tHen find the 
number of levels which would yield the largest power for these values. 
In each case the possible number of levels ranged from two to N/2 
since at L- 1, can not be estimated, and at L = N neither 
nor a|'can be estimated. 

For example, in constructing Table 2 where 6^ = 0.00, 9= 0.60, 
and N = 20 the following values wer^ generated for the levels from two 
through ten (N/2): Levels (L) ppwer 

2 .43750 

3 .46739 
* ,50688 
5 .48910 



The author is indebted to Miss Ruth A. Wagenhofer for her effort in 
compiling the tables. 
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Table 1 









Power Values and Optimum Number 












of 


Levels for Rixed Total Numbers 




• 








of 


Observations in the One-Way 












Random-Effects Analysi 


s of Variance 
















a = .01 




















= 0.00 










e/N 




0.2 


0.4 


0.6 


0.8 


1.0 


2.0 


3.0 


4.0 


10 


2. 


.045 


2,. 089 


2, . 132 


2 M7 

^ , . X / 4^ 


2 2nR 


2 1^1 
^ , . J4X 


9 496 


2 AftS 


20 


2, 


.114 


2 ?14 




L '^OA 

4 J . J 


4 , • 4 / X 


7nA 

J, . /UO 


5, .706 


5,. 881 


30 


3, 


.180 


3 .348 


3 348 


^ , . JOO 


J , « UUO 




6,. 875 


10, .977 


40 


3 


- 246 




4 463 


S 71 
J f i Xu 


ft 7QS 


xu , . 7 J4 






50 


3, 


.306 


5 561 


S S61 


7 AOQ 


in ft7ft 


xu, . 7ox 


10,. 981 


16,. 999 


60 


4, 


.383 


6 666 


6 644 


1 0 ftfti 
XU , . oox 


xu , . 7 JU 


Xj , . 77'» 




19 1 nn 
xz, X. uu 


80 


5, 


.439 


7 713 


7 71 3 

i f • t XJ 


1 0 Q9A 


X4 , • 70U 


1 A QQQ 
X** , . 770 


14,. 998 


10,1.00 


80 


5, 


.497 


8 770 


8 770 


10 Q^O 


Ifi 077 
Xu , . 7 / / 


7n QQQ 
ZU, . 777 


20,. 999 


9,1.00 


90 


5, 


.547 


Q 817 


Q 81 7 

7 9 « OX/ 


1 C Q7n 
XD , . 7 /U 


1 Q QQQ 
XO , . 700 


X J , X. uu 


15,1.00 


8.1.00 


100 


5, 


.591 


XV/ 9 • O ^ ^ 


10 8SS 


90 QftO 


2n QQT 


X4,x. uu 


14,1.00 


8,1.00 


300 


15, 


.965 


30,. 999 


30, .999. 


7,1.00 


6,1.00 


6,1.00 


6,1.00 


5,1.00 


500 


15, 


1.00 


8,1.00 


8.1.00 


7,1.00 


6,1.00 


6,1.00 


6,1.00 


5,1.00 












= Q.IO 










e/N 


0.3 


0.5 


0.7 


0.9 


1.0 


2.0 


3.0 


4.0 


10 


2, 


.032 


2, .059 


2, .039 


2,. 118 


2, .132 


2, .250 


2, .334 


2, .396 


20 


2, 


.057 


4,. 120 


4,. 194 


4, .267 


4, .302 


.5, .567 


5, .718 


5,. 804 


30 


3, 


.082 


5,. 188 


5, .302 


6,. 405 


6, .454 


6, .751 


10,. 891 


10, .948 


40 


4, 


.106 


5,. 252 


8, .395 


8, .527 


8, .582 


10,. 878 


10,. 958 


13, .984 


50 


5, 


.130 


7,. 309 


10,. 484 


10,. 629 


10,. 686 


10., 931 


16,. 984 


16,. 996 


60 


6, 


.153 


10,. 371 


10,. 569 


12,. 613 


12,. 768 


15, .971 


20,. 996 


20,. 999 


70 




.176 


10,. 430 


14,. 636 


14,. 781 


14,. 830 


23,. 985 


23,. 999 


17,1.00 


80 


8, 


.199 


10,. 479 


16,. 697 


16,. 832 


16,. 877 


20,. 994 


26,1.00 


15,1.00 


90 


10, 


.223 


15,. 530 


15,. 751 


18,. 876 


18,. 912 


30,. 997 


18,1.00 


14,1.00 


100 


10, 


.245 


14,. 569 


20,. /95 


20, .907 


20, .978 


25, .999 


16,1.00 


12,1.00 


300 


30, 


.632 


50,. 967 


60,. 998 


50,1.00 


30,1.00 


17,1.00 


10,1.00 


8,1.00 


500 


62, 


.850 


50,. 998 


22,1.00 


16,1.00 


15,1.00 


10,1.00 


8,1.00 


8,1.00 



ERIC 



Table 1 continued 



-8 



e 0.50 



e 8.7 



N 

10 

20 

30 

40 



0.8 



2, 
4, 
5, 
8, 



50 10,. 

60 12, . 

70 14,. 

80 16,. 

90 18,. 
100 20, . 
500 125,. 
1000 250,. 



018 2,. 023 
024 4,. 034 
029 6,. 043 
034 8,. 053 
039 10,. 062 
043 12,. 071 
048 14,. 081 
052 16,. 090 
056 18,. 099 
06L 20,. 109 
255 125,. 511 
503 250,. 828 



0.9 

2, .028 
4, .045 
6,. 060 

10,. 091 
ia,.107 
14,. 122 
20,. 139 
18,. 154 
25,. 175 
125,. 739 
333, .966 



1.0 . 

2,. 034 

4,. 057 
6,. 080 
8,. 102 
10,. 125 
15,. 149 
14,. 171 
20, .197 
18,. 218 
25,. 246 
125,. 884 
190,1.00 



!.0 



2,. 095 

5,. 218 
10,. 329 
10,. 445 
16, .526 
20,. 634 
23,. 703 
20,. 762 
30,. 823. 
33^.^856 
68,1.00 
35,1.00 



3.0 

2, .155 

5, .380 
10,. 570 
13,. 344 
16',. 799 
20,. 885 
23,. 926 
26, .953 
30, .975 
33,. 985 
20,1.00 
20,1.00 



«.0 

2,. 208 

5,. 508 
10,. 730 
13,. 847 
16,. 916 
20,. 964 
23,. 982 
26,. 991 
30,. 996 
33,. 998 
15,1.00 
15,1.00 



5.0 

3, .259 

5, .604 
10,. 828 
13,. 919 
16, .963 
20,. 988 
23, .995 
40, .998 
45, .999 
30,1.00 
13,1.00 
13,1.00 



e 1.2 



N 
IC 
20 
30 



2,. 015 
4,. 017 
6,. 020 
40 10,. 022 
50 10,. 024 
60 15,. 026 
70 14,. 027 
80 20,. 029 
90 22,. 031 
100 25,. 033 
500 166,. 103 
1000 333,. 202 



1.4 

2,. 020 
5,. 027 
6,. 034 
10,. 040 
10,. 046 
15, .053 
17, .058 
20, .065 
30, .071 
25, .078 
166,. 365 
333, .677 



1.6 

2, .035 
5,. 039 
6,. 032 
10,. 065 
10,. 07/5 
15,. 091 
23,. 102 
20.. 118 
30,. 132 
25,. 145 
166,. 678 
333,. 944 



e ^ = l.CO 
1.8 



2, .032 
5,. 054 
6, .073 
10,. 096 
12,. 113 
15,. 140 
23,. 161 
20,. 185 
30,. 211 
33, .233 
166,. 882 
250,1.00 



2.0 


3.0 


4.0 


5.0 


2, .038 


2, 


.074 


2, .111 


2,. 146 


55.070 


5, 


.166 


5, .270 


5,. 365 


10,. 098 


10, 


.260 


10, .429 


10,. 570 


10,. 132 


10, 


.344 


13,. 550 


13,. 702 


12, .157 


16, 


.425 


16, .653 


16,. 799 


20,. 198 


20, 


.525 


20,. 760 


20,. 885 


23,. 229 


23, 


.592 


23, .821 


23,. 926 


20,. 262 


26, 


.651 


26,. 869 


40, .955 


30,. 302 


30, 


.721 


30,. 914 


30,. 975 


33, .333 


33, 


.765 


33,. 938 


50,. 985 


166,. 967 


94, 


1.00 


34,1.00 


24,1.00 


142,1.00 


52, 


1.00 


34,1.00 


24,1.00 











Table 2 
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Power Values and Optimum Number 










of 


Levels for Fixed Total Numbers 










of 


Observations in the One-Way 








• 


Random-Effects Analysis of Variance 














0 « .05 


















0Q= 0.00 












0.2. . 


0.4 


0.6 


0.8 


1.0 


2.0 


3.0 


4.0 


10 


2,. 142 


2,. 220 


2,. 282 


2,. 333 


2,. 374 


3,. 518 


3, .622 


5, .693 


20 


2,. -241 


4,. 386 


4, .507 


4, .596 


4, .662 


5, .847 


5, .914 


5, .945 


30 


3,. -342 


5,.53P 


5, .666 


§..756 


,6, .818 


10, .949 


10,. 984 


10, .994 


40 


4 .424 


5,. 645 


'8, .768 


8,. 854 


8,. 904 


10, .984 


1 1 Odd 


ij , . yyy 


50 


5, .495 


5, .724 


10,. 843 


10,. 914 


10,. 950 


16, .994 


16, .999 


12,1.00 


60 


5, .564 


6, .792 


10,. 900 


12,. 950 


12,. '974 


20, .999 


12,1.00 


10,1.00 


70 


5 f>7^ 


7, .843 


10, .934 


14, .171 


14,. 987 


23,1.00 






80 


5, .668 


10, . 883 


10,. 955 


16, ,983 


16, .993 


13,1.00 


9,1.00 


8,1.00 


90 


6 .71 ^ 


10,. 913 


15, .965 


18,. 997 


18,. 997 


11,1.00 


o, J . 00 


/,1.00 


100 


7, .746 


10,. 933 


14,. 980 


20,. 995 


20,. 998 


10,1.00 


8,1.00 


7,1.00 


300 


20,. 989 


23,1.00 


8,1.00 


7,1.00 


6„1.00 


5,1.00 


5,1.00 


5,1.00 


500 


10,1.00 


7,1.00 


7,1.00 


6,1.00 


6,1/00 


5,1.00 


5,1.00 


5,1.00 










00 = 0.10 










. 


0.3 


0.5 


0.7 


0.9 


1.0 


2.0 


3.0 


4.0 


10 


2,. 112 


2,. 169 


2,. 220 


2,. 263 


2, .282 


3, .436 


3, .547 


5, .628 


20 


4,. 163 


5,. 283 


4, .386 


5, .473 


5,. 512 


5, .753 


5,. 854 


10,. 907 


30 


5, .211 


5, .377 


6, .513 


6,. 618 


6,. 661 


l-,.884 


10,. 962 


10,. 984 


40 


5,. 255 


8,. 445 


8, .616 


8,. 727 


10, .771 


10, .952 


13, .988 


13, .996 


50 


7,. 290 


10,. 526 


10, .698 


10,. 806 


10,. 843 


16, .978 


16,. 996 


16, .999 


60 


10,. 325 


10,. 592 


12, .764 


15,. 863 


15,. 897 


20, .993 


20,. 999 


14,1.00 


70 


10,. 364 


10,. 644 


14, .816 


14,. 904 


14, .930 


23,. 997 


17,1.00 


13,1.00 


80 


10,. 397 


16, .692 


16,. 857 


20,. 934 


20,. 955 


26,. 999 


15,1.00 


11,1.00 


90 


10,. 426 


15, .638 


18, .890 


18,. 954 


13,. 969 


30,1.00 


14,1.00 


10,1.00 


100 


11,. 453 


20,. 771 


20,. 915 


25,. 969 


25, .981 


20,1.00 


12,1.00 


10,1.00 


300 


37, e 821 


50, .991 


48,1.00 


26,1.00 


21,1.00 


15,1.00 


8,1.00 


7,1.00 


500 


71,. 949 


29,1.00 


17,1.00 


13,1.00 


12,1.00 


9,1.00 


8,1.00 


7,1.00 



ERIC 



Table 2 continued 



10 



N 






10 


2, .076 


2, .090 


Ort 
zU 


5, .095 


5, .121 


on 


6, .109 


6, .144 


/.n 


8, .122 


10,. 167 


cn 
jO 


IO9.I34 


10,. 186 


An 


1 c 1 / c 

15^.145 


15, .207 


7n 
/U 


1 A ICC 

!4,.155 


14, .224 


fin 


^U, .166 


20,. 245 


on 


lo, .175 


18,. 260 


1 nn 


23,. 186 


25,. 28 3. 




123). 499 


125,. 749 




250,. 745 


333,. 944 




e 1.2 


1.4 


N 






10 


2, .065 


2,. 081 


20 


5, .076 


5,. 105 


30 


10,. 083 


10,. 123 


40 


' 10,. 090 


10, .139 


50 


12, .094 


16, .152 


60 


15, .101 


20,. 170 


70 


23,. 106 


23,. 183 


80 


20, .112 


20,. 196 


90 


30,. 116 


30, .212 


100 


25, .121 


33,. 224 


500 


166,. 277 


166, .625 


1000 


333,. 438 


333,. 867 



8^ - 0.50 

0 



0.9 1.0 2.0 

2, .103 2, .116 3, .239 

5, .147 5, .175 5, .429 

6, .181 6, .218 10, ,578 

10,. 215 10, .265 13, .677 

10,. 242 10, .300 16, .755 

15, .275 15, .344 20, .832 

14,. 298 17, .373 23, .875 

20, .330 20, .415 26, .907 

22,. 351 30,. ■•2 30, .938 

25,. 381 J5,.480 33, .955 

125, .899 166, .967 27,1.00 

237,1.00 115,1.00 27,1.00 



e =1.00 

0 



1 


.6 


1.8 


2.0 


2, 


.096 


3, .112 


3, .129 


5, 


.137 


5,. 169 


5, .203 


10, 


.168 


10, .216 


10,. 267 


10, 


.195 


10,. 255 


13, .317 


16, 


.219 


16, .291 


16,. 365 


20, 


.251 


20,. 338 


20,. 424 


23, 


.274 


23, .371 


23,. 466 


26, 


.296 


26,. 402 


26,. 505 


30, 


.325 


30,. 442 


30,. 533 


33, 


.346 


33, .471 


33,. 587 


165, 


.868 


166,. 967 


166,. 993 


333, 


.992 


142,1.00 


100,1.00 



3.0 4.0 5.0 

3, .343 5, .429 5, .509 

5,. 601 10, .720 10,. 809 

10,. 784 10,. 884 10,. 934 

13i.870 20,. 945 20,. 978 

16, .924 25, "7/ ?5,.993 

20,. 963 20.. 9)1 i0..998 

23, .979 35, .996 j5,.999 

26, .988 40, .999 24,1.00 

30, .994 45, .999 21,1.00 

33, .997 jO,1.00 20,1.00 

16,1.00 13,1.00 11,1.00 

16,1.00 13,1.00 11,1.00 



3.0 4.0 5.0 

=3,. 209 3, .281 5,. 350 

5,. 361 10, .491 10,. 610 

10, .501 10, .672 10, .784 

13, .595 13, .771 20, .878 

16, .672 25, .848 25, .935 

20, .756 20, .906 30, .966 

23,. 805 35,. 938 35,. 982 

26, .845 40, .961 40,. 991 

30,. 887 45,. 976 45,. 996 

33, .911 50, .985 50, .998 

77,1.00 25,1.00 19,1.00 

41,1.00 25,1.00 19,1.00 
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Table 3 

Power Values and Optimum Number 
of Levels for Fixed Total Numbers 
of Observations in the One-Way 
Random- Effects Analysis of Variance 
a » .10 



e » 0.00 

0 





e 0.2 


0.4. 


0.6 - 


0.8 


. 1.0 


, 2.0 


'3.0 


4.0 


N 
10 


2,. 225 


2, .314 


2,. 380 


2,. 243 


2, .470 


3,. 633 


5,. 737 


5,. 808 


20 


4, .331 


4, .501 


4, .615 


5,. 694 


5, .755 


5, .897 


5, .944 


10,. 972 


30 


3,. 444 


5, .637 


6, .754 


6, .830 


6, .877 


10, .973 


10,. 992 


10,. 997 


40 


4, .531 


5, .734 


P,.843 


8, .906 


10, .941 


10, .992 


13,. 998 


13,1.00 


50 


5, .602 


7, .801 


10,. 900 


10,. 948 


10, .971 


16, .997 


16,1.00 


10,1.00 


60 


5, .662 


10, .855 


10,. 938 


12, .971 


15, .986 


20,. 999 


10,1.00 


8,1.00 


70 


7, .711 


10, .897 


10,. 960 


14, .984 


14, .993 


14,1.00 


9,1.00 


7,1.00 


80 


8, .753 


10, .926 


10, .975 


16, .991 


20, .997 


10,1.00 


8,1.00 


7,1.00 


90 


6, .791 


10, .945 


15,.98'v 


18, .995 


18, .998 


9,1.00 


8,1.00 


6,1.00 


100 


9, .821 


11, .959 


20,. 990 


20,. 998 


20,. 999 


9,1.00 


7,1.00 


6,1.00 


300 


20, .994 


20,1.00 


7,1=00 


7,1.00 


6,1.00 


5,1.00 


5,1.00 


5,1.00 


500 


9,1.00 


7,1.00 


6,1.00 


6,1.00 


6,1.00 


5,1.00 


5,1.00 


5,1.00 





6 0.3 


0.5 


N 






10 


2, .188 


2, .258 


20 


4, .258 


4, .396 


30 


5, .316 


6, .496 


40 


5, .363 


5, .558 


50 


7,. 406 


10,. 646 


60 


10,. 448 


12, .702 


70 


10, .487 


14, .750 


80 


10,. 520 


16, .790 


90 


15,. 550 


15, .824 


100 


14, .577 


20,. 853 


300 


50,. 892 


60, .996 


500 


71, .975 


23,1.00 



9 = 0.10 

0 

0.7 0.9 1.0 

■ 2, .314 3, .360 3, .385 

5, .505 5, .592 5, .628 

6, .628 6, .720 6, .754 

8, .721 10, .816 10, .850 

10, .791 10, .873 10, .900 

12, .844 15, .918 15, .940 

14, .883 14, .943 14, .960 

16, .913 20, .964 20, .977 

18, .935 18, .975 22, .984 

20, .952 25, .984 25, .991 

33,1.00 23,1.00 21,1.00 

15,1.00 12,1.00 11,1.00 



2.0 


3.0 


4.0 


3, .558 


5, .678 


5,. 758 


5, .828 


10,. 907 


10,. 953 


10,. 939 


10,. 980 


10,. 992 


13, .974 


13,. 994 


20,. 998 


16, .989 


16, .998 


25,1.00 


20,. 997 


20,1.00 


12,1.00 


23, .999 


14,1.00 


10,1.00 


26,1.00 


13,1.00 


10,1.00 


18,1.00 


11,1.00 


9,1.00 


16,1.00 


11,1.00 


9,1.00 


10,1.00 


8,1.00 


7,1.00 


8,1.00 


7,1.00 


7,1.00 



Table 3 continued 
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e 


0.7 


0.8 


N 






10 


2,. 140 


2,. 158 


20 


5,. 169 


5, . 205 


30 


b,.189 


6,. 236 


40 


10,. 208 


10,. 268 


50 


10,. 223 


10, . 291 


60 


15,. 240 


15,. 319 


70 


14,. 253 


17,. 338 


80 


20, . 268 


20,.J65 


90 


18,. 279 


22,. 382 


100 


25,. 294 


25,. 407 


500 


125,. 636 125,. 845 


1000 


250,. 840 


333,. 976 



= 0.50 

0.9 1.0 2.0 

3, .177 3, .196 3, .355 

5,. 241 5,. 276 5,. 551 

10,. 284 10,. 333 10,. 700 

10,. 328 10,. 386 13,. 784 

10,. 358 12,. 423 16,. 845 

15,. 399 15,. 474 20,. 901 

17,. 424 23,. 507 23,. 930 

20,. 460 20,. 548 25,. 950 

30,. 485 30,. 580 30,. 969 

25,. 515 25,. 612 33,. 978 

125,. 947 166,. 985 23,1.00 

142,1.00 93,1.00 23,1.00 



3.0 


4.0 




5.0 


5,. 483 


5, .582 


5 


,.657 


10,. 720 


10,. 829 


10 


,.892 


10,. 863 


15,. 933 


15 


,.968 


13, .925 


20,. 974 


20 


,.991 


25,. 961 


25, .990 


25 


,.997 


20,. 982 


30,. 997 


30 


,.999 


23,. 990 


35, .999 


23 


,1.00 


40,. 995 


40,1.00 


20 


,1.00 


30,. 998 


30,1.00 


18 


,1.00 


50,. 999 


25,1.00 


20 


,1.00 


14,1.00 


12,1.00 


10 


,1.00 


14,1.00 


12,1.00 


10 


,1.00 



e 


1.2 


1.4 


M 








10 


3, 


.124 


3,. 1.49 


20 


5, 


.141 


5,. 184 


30 


10, 


.152 


10,. 211 


40 


10, 


.162 


10,. 232 


50 


15, 


.171 


16,. 253 


60 


20, 


.180 


20,. 275 


70 


23, 


.187 


23,, 292 


80 


20, 


.194 


26,. 308 


90 


30, 


.202 


30,. 328 


100 


33, 


.209 


33,. 343 


500 


166', 


.408 


166,. 759 



1000 333,. 580 333,. 928 



Qq = 1.00 

1.6 1.8 2.0 

3, .173 3, .196 3,. 219 

5, .227 5, .269 5, .310 

10, .272 10, .333 10, .392 

13,. 305 13,. 379 13,. 449 

16,. 340 16,. 426 16,. 506 

20,. 375 20,. 472 20,. 561 

23,. 402 23,. 507 23,. 602 

26,. 427 26,. 540 26,. 639 

30,. 459 30,. 580 30,. 683 

33,. 482 33,. 608 33,. 613 

166,. 929 166,. 985 166,. 997 

233,1.00 116,1.00 82,1.00 



3.0 4.0 5.0 

3, .330 5, .422 5, .502 

10, .492 10, .639 10,. 743 

10, .632 10, .779 15, .871 

13,. 716 20,. 864 20,. 936 

16,. 787 25,. 918 25,. 969 

20,. 846 30,. 951 30,. 985 

23, .883 35, .971 35, .993 

26,. 913 40,. 983 40,. 997 

30, .938 45, .990 45, .999 

50,. 953 50,. 994 50,. 999 

69,1.00 22,1.00 17,1.00 

36,1.00 22,1.00 17,1.00 
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Levels (I) Power 

6 .40543 

7 .23867 

8 .26410 

9 .28892 
10 .31314 

Maxinum power (.50688) is reached when there are four levels with 
five (N/L) elements under each level. The values entered in the 
table were 4, .507. 

Tables 1, 2, and 3 have alpha levels of .01, .05, and .10 
respectively. Each of the tables has set at 0.00, 1.10, 0.50, 
and 1.00, with selected values of and N. These values were • 
selected so as to provide information for common one-way random-effects 
designs. 

The power values in the tables have been rounded to three decimal 
places. In certain situations sample size was so large that power 
exceeded .9995 for many values of L. When this happened, the first 
value of L whose corresponding power exceeded .9995 was selected as 
being the optimum number of levels to use in the analysis of variance. 
Since optimum L generally increases with Increasing values of 9 or N 
the preceding process explains why some of the values of L decreased with 
increasing 9 or N. 

One should be aware of the fact that the tables presented here 
provide "optimal" sample size and number of levels. Given specified 
parameters, one can not take different values of L and find larger 
power. However, an e:<periment having slightly less power might be 
satisfactory if the number of levels is limited by cost or some other 



factor. For example, if a = .10, » 0.10, N » 60 and 6 = 1.0, 
maximum power (.940) is found when L » 15 and n = 4. Consider the same 
parameters -k^en less than maximum power is desired, then with 
L at 10, n = 6, power »,.928; with L = 5, n = 12, po^er = .845. Such 
power values may be found for situations not covereJ in the tables 
presented by using extensive F tables such as those found in 
Graybill (1961) or by using a computer program described in - 
Barcikowskl (1972) , the mainline of which is in appendix A. 

Example 

The following discussion is based on the problem mentioned 
earlier concerned with teaching first graders how to read. Four 
different cases will be presented with the availability of teachers 
(levels) and students (observations) varying in each case. 

Suppose that previous research has indicated that teachers 
under conditions similar to those in this experiment are about 
one-tenth as variable as their students; i.e., 6^ = 0.10. Based on 
the investigator's experience, e.g., a pilot study, he has decided 
that it would not be worth the effort to use the method of teaching 
reading if the ratio of the variability among teachers to the 
variability among students is 1.0 or larger. This procedure leads 
him to the following hypothesis and its alternative: 

"o= 4 1 -lOOe 
■■■■ 2 2 
Hq: oa > -lOOe 

That is, the reading method will not be adopted for use if H 

0 ' 

is rejected, but will be adopted if the researcher fails to reject 

Hq. 
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Power for this experiment is based on the fact that if 
e >^ 1.0 we would like to detect this with high probability. Armed 
with this information the researcher may now turn to Tables 1, 2, 
and 3 and select the number of teachers, L. and the total number 
of students, N, which will keep the probability of making a Type 
II error small. His selection of an o level may be dependent on 
the availability of teachers and/or students to achieve a certain 
power. However, let us assume he has selected a = .05, and then 
consider the following situations which will influence his selection 
of a power value. 

Case I; L and N Plentiful 

Optimally the researcher is in a position where he may select 
anf number of students and teachers. If this is the case he would 
turn to Table 2 under 6^ = .10, 8 = 1.0 and find the power desired 
for this. experiment. If he selectao power at .98 then he might use 
25 teachers with each teacher teaching 5 pupils. However, if he 
selected power at 1.00, 12 teachers with 41 pupils 500/12 each might 
also be selected. Here one would not need to use N = 500 pupils in 
order to achieve maximum power but only 12 X 41 * 492. 
Case II; L and N Fixed 

Consider the problem where due to costs there are a fixed 
number of teachers and pupils. Let the maximum number of teachers 
and students that can be sampled at random from their respective 
populations be 15 teachers an^O pupils. Then Table 2 indicates 
that using only 10 teachers ^■5 pupils each will yield an experiment 
having maximum power of .843flPere using the .extra teachers will 
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only result In a statistical test having less power.' 

If there were a limit of 20 fiachers and 75 students then the 
power of the statistical test using the optimum number of levels- 
would lie between .930 and .955. Here 15 teachers would probably 
be a good choice since 15 divides into 75 evenly and is in the 
range of number of levels, 14-20, considered optimal for 6 = 1.0. 
However, one would certainly want to check this by hand calculation. 
Case III; L Fixed; n Plentiful 

If there are a fixed number of teachers and plenty of students 
the researcher may simply look in Table 2 at Qq =0.10 under 
e = 1.0 until power appeared to be sufficient. If the investigator 
desired power at . 80 and he had 15 teachers then he could use 10 
teachers and 5 students per teacher. If he had fewer than 3 teachers 
and desired power at .80 then he would have to do his calculations 
by hand. 

Case IV; L Plentiful; M Fixed 

If the number of students is fixed at 80 then the optimum 
number of teachers to select would be 20. In this case. the power 
of the experiment would be .955. Here again, the researcher may 
not be satisfied with the student-teacher ratio, but .955 is the 
largest power he can achieve using 80 students. If he is willing 
to have less power, then using hand caculations he will find using 
8 teachers with 10 students each results in a power of .906 or 5 
teachers with 16 students each results in a power of .'821. 
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The Two-Way Design 

In the two-way random-effects deoign, one may test three 
possible null hypotheses. They are the hypotheses concerned 
with the variance of the main effects, factors A and B, and 
the variance of the interaction effects between factors A and 
B. The interaction hypothesis would only be tested provided the 
number of observations in a cell was two or more. In this presentation 
it is assumed that the. number of observations in eaclj cell are equal. 

The null hypothesis concerned vil:h factor A would be written 

0A 1 ^oA 

with the alternative 

The null hypothesis concerned with factor B would be written 
with the alternative 



%2* ~ ^OB' 



ERIC 



The null hypothesis concerned with the interaction of factors 
A and B would be written 

Where 9^ = aj/ (ag-fna^g) ; Gg = a^/(a^4<T^); e^^ = alg/a^, and 

^oA» ^0B» ^o^j^ are preassigned constants, greater than or equal 

to 2ero. Here the variables are defined as follows: 
2 

= the variance of the effects for factor A, 

2 

ag - the variance of the effects for factor B, 

2 

^AB variance of the interaction effects, 

2 

Oe = the error variance, 
n e the number of observations in a cell of the design. 
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The power of the F test (Scheffe, 1954) on K , , ,H „ and H may be 

OA oB OAB ' 

calculated using 

Power = PrCFd-LV^g) > •F(cf;I-l,V^g) (l+Jne^A)/(H JneA)}; 

Power = Pr{F(J-l,V^g) 1 F(a;I-l,V^g) (l+Ine^g)/a+Inf g)} 

Power = Pr{F(V^g.V^) F(«;V^.V^) (l+ne^^)/(l+v.e^.,>. 
respectively. Here the parameters are defined as: 

F = the F statistic, 

a = the level of significance, 

I *= the niimber of levels in factor A, . 

J = the niimber of levels in factor B, 

V = ij(n.l)r 
e 

V^-= (I-1)(J-1) 

Using the preceding power formulas with the central F distribution 

a computer program^ '(Appendix B) was written to calculate power for the 

statistical F tests used to test the three possible null hypotheses • The 

procedure used was to select a, 6, e^and the total number of observations 

(N) and then find the number of levels which would yield the largest 

power for these values across the tests of H , H „, and H . In each 

oA oB oAB 

case the number of levels of B (or A) ranged from 2 to N/4 while the 
number of levels of A (or B) ranged from N/4 to 2. This procedure was 

used since if I or J was one, then either o. or c can not be estimated. 

A • B ^ 

and if I or J was greater than 1^/4, 4 can not be estimated. 

AB 

^One minor limitation of this program is that it assumes 
®A " ®B " %B ^OA " ^oB " ^oAB- ^^^^^ different values of 

06 andce^ one can bypass this limitation by submitting the program 
once for each of the different pairs of G and 6 , and then 

.making comparisons of levels across programs. * 
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Example 

Tables 4, 5 and 6 provide an example of the optimum number 
of levels to test 9^; 9^ and 6^. respectively. In all three 
tables a = .05. 9 = .80. 9^ = 0.0 and N = 40. The selection of 
the preceding parameters would prevent the investigator with a 
problem. In table 4, power is optimal .777. for testing H 
when there are two observations in a cell, I = 5 and J = 4. in 
table 5. power is optimal. .777. for testing H ^ when there are 
two observations in a cell. I = 4 and J = 5. That is. maximum 
power is achieved for testing on each of the factors when the 
number of levels is five. This Is not a major problem since in 
both cases the power for testing H or slips to .764. So 
that, if one tested H^^ with I = 4 and J = 5 power would be .777 
for this test and .764 for the test of K^^. m table 6 the test 
°^ "OAB ^^^^^^^s maximum power. .638. when there are four ob- 
servations in a cell and J = 2. I = 5 or I = 2 ..ad J = 5. Under 
these conditions the power of the test of main effects drops to 
.538 and .555 depending the number of levels selected for factors 
A and B. 

In this example it would seem that one would be forced to in- 
crease sample size and then test the interaction hypothesis. If 
one failed to reject the interaction hypothesis and if power was not 
sufficient for the main effects test, a second experiment would 
have to be run with d^lfferent levela to test the main effects. 
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Although in practice one would probably have different values 
°f ^A' ^"'J %y with different 's, the preceding example 
provides some indication of the dil Ticulties one might encounter 
in selecting optimum numbers of levels to arrive at optimum power 
..1 the two-way random-effects analysis of variance. 

The Nested Design 

In considering nested designs these are as many hypotheses to 
be tested as there are factors. In this design interactions are 
assumed to be nonexistant. The design discussed here is composed 
of three factors, A, S and C, from which the reader may generalize 
to other nested designs. There are an equal number of observations, 
n, in each cell and the number of levels nested within a level of 
a factor are equal. There are I levels of A; J levels of B nested 
within each level of A; K levels of C are nested within each level 
of B. Figure L provides an example where I = 2, J = 3, k = 12, 
and n = 4. In figure 1 there are A subjects in each cell, 3 levels 
of B are nested within each level of A and two levels of C are 
nested within each level of B. 



Figure 1 
An Example Nested Design 
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The null hypothesis concerned with each factor would b« 
written 



H : e < e 
ox X - ox 



with the alternative being 

H. : < e 
Ax X — "ox 

where x may take on the factor letter designations A. B, or C. 
Where 

2a 2 
A » " - 



2 

a 

e 



^ox Pi^eassi^ned constants. 

The power of the F test (Scheffe, 1954) on H . H .and 

^oC calculated using 

Power » Pr{F(I.l,V ) > F(o;I«l,V )(l+JKn0 )/(l+JKn0 )}; 

^ oA A 

Power » Pr{F(Vg.B^) ^F(a;Vg.V^) (1+Kne^^)/(1+Kneg)}; 

Power = Pr{F(V^.VB) >Ha;V^,V^) (l-^e^^) / (l+nQ^)} , 

respectively. Here the parameters are defined as: 

F = the F statistic, 

a = the level of significance, 

I « the number of levels in factor A, 

J = the number of levels in factor B, nested within each 
level of A, 

K - the number of levels in factor C, nested within each 
level of B, 

Vb=I(J-1). V -JI(K-l), V » IJK(n-l). 
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Given a fixed and Q^^, the same problems encountered In the 
two-way analyses will be encountered with nested designs. That 
Is, the power for each factor varies and is partially dependent 
on the number of levels In the other factors. Here the problem 
Is even more complex since there Is a large number of designs 
which can be generated given flexibility in sample size and 
number of levels for each factor. 

Summary 

This presentation furnishes researchers in the social sciences 
with procedures for selecting sample size and optimum number of 
levels for the random-effects analysis of variance. Hopefully, in 
so doing, it has brought to the attention of these resaarchers 
problems which could be investigated using this analysis of variance 
model. 



26 



REFEllEHCES 



^'""'"if-'S"^' ^* ^ computer program for selecting optimum 

^°^'"ccic/*r''°''^%'^ "^Pl^ determination for 

aTrs" ^ho'LT'r m"' ^"°-^Plf hypothesis tests. Educational 
■2 i£. J^s. 'etiological Measurement . 1970, 30, 811-831. 

^°^'"n4*Yo1-^T^^S^.^ analjrsis for the behavioral sciences . 
I w loi.k: Academic Press, 1969. - 

^^''""anLs;: * v^'T* J' ^" Introduction to statistical 
anal2sj£. New York: McGraw-Hill Book CompaT^,T969: 

^^^''and*Dl;:hol^'^"^'^' ^' Statistical me^hoH■. -r. .^w^^ 

"^^^'RLL^^4^^-^£5Z£holo.Si^ New York: Holt, 
Rinehart, and Winston, 1965. ' 

^'"■scie^ie.^^lfSr^r^m^. £rocadurss for the beMvloral 
i^sS:"l968^ • '^"1"""^== Brooks/Col. PubliSHSi 



27 



Appendix A 



•Only the MAINLINE of the progran is included here. The subroutines 
FSTAT and PLEVEL are very long and were excluded. FSTAT returns 
the F statistic (F) given alpha (ALPHA) and degrees of freedom 
(Dl, D2); PLEVEL returns the probability of a given F statistic 
(P) given degrees of freedom (Dl, D2). Similar routines should 
be available at most computer centers. 
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SSuc^'T" EJiT.-CTS-POWER' WAS PR0G8A^^MED BY ROBERT S. 
BARCIKOWSKI AT OHIO UNIVERSITY 1971-1972 

SSsSS/f^'^^f'a'^^^lJ!-,,^"^^"^^^ Al^n) PSYCHOLOGICAL 
PROBLEM CARD 1 

COLUMNS 1-4 NIHBER OF DATA SETS TO BE READ (I.E.. THE NUMBER OF 
TIMES PROBLEM CARDS 2-4WILL BE REPEATED) . 
(14 FORMAT) 



PROBLEM CARD 2 



COLUMNS 1-3 

4-6 
7-9 
10-43 

PROBLEM CARD 3 



ALPHA, LEVEL OF SIGNIFICANCE FOR THE F TEST. 
(F3.2 FORMAT) 

NUMBER OF PHI'S TO BE READ. (13 FORMAT) 
NUMBER OF N'S TO BE READ. (13 FORMAT) 
PHI (SUB 0), HYPOl-HESIS VALUE. (F4.2 FORMAT; 



COLUMNS 1-3 
4-6 



FIRST VALUE OF PHI. 

SECOND VALUE OF PHI, ETC., UP TO 20 VALUES. 
(20F3.1 FORMAT) 



PROBLEM CARD 4 



COLUMNS 1-4 
5-8 



FIRST VALUE OF N. 
SECOITO VALUE OF N, 
(20F4.0 FORMAT) 



ETC., UP TO 20 VALU2S. 



RE-POCOl 
RE-P0002 
RE-P0003 
RE-P0004 
RE-P0005 
RE-P0006 
RE-P0007 
RE-P0008 
RE-P0009 
KE-POOlO 
RE-POOll 
RE-P0012 
RE-P0013 
RE-P0014 
EE-P0015 
RE-P0016 
RE-P0017 
RE-P0018 
RE-P0019 
RE-P0020 
RE-P0021 
RE-P0022 
RE-POO 23 
PvE-P0024 
RE-P0025 
RE-P0026 
RE-P0O27 
RE-P0028 
RE-P0029 
RE-P0030 
RE-P0031 
RE-P0032 
RE-P0033 



K^ut^.u FORMAT) 'vwi 

DOUBLE PRECISION D1,D2,ALPHA.F FPOW P RE-P0036 
Bt7An ' » "«,r.rruw,r RE-P0037 
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READ 11, NJET 
F0RMAT(I4) 
DO 10 JET - 1,NJET 
READ 7,ALPHA,NPHIS,NrOT,PHO 
F0RMAT(F3.2,2I3,F4.2) 
READ 1,(PHI(J),J-1,NPHIS) 
FORMAT (20F3.1) 
READ 6,(T0TNUM(J),J-1.MT0T) 
FORMAT (20F4.0 
lACC - 5 

PRINT 8,NT0T,NPHIS,ALPHA 



RE-P0038 
RE-P0039 
RE-P0040 
RE-P0041 
RE-P0042 
RE-P0043 
RE-P0044 
RE-P0045 
R3-P0046 
RE-P0047 
RE-P0C48 



CLUES OF PHI, ir-FsM/' '1^' " VARE-POOSO 

RE-P0052 
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DO 4 T.l - 1,NT0T RE-P0053 
NU-'ITCr^TOTNUMdl) RE-P0054 
fRiW. 9,NUyiT0T,PH0 RE-P0055 
FOR'U-.TdHO' THE TOTAL m-IBER OF OBSERVATIONS IS 'l6' THE NULL HYP0RE-P0056 
CTHE31S HAS A VALUE AT 'F5.2) RE-P0057 
JE?]i><lOTI]UM(Il)/2.0 RE-P0058 
0.0 A..I2 = l.NPKIS RE-P0059 
PRIRT 5, PHI (12) RE-P0060 
FOa^'^iXlKO' ************************************************** '^^ RE-P0061 
C Zm PHI COEFFECIENT IS 'Fe.a,//' ***************ft***************RE-P0062 
C***'r.. A;****.********* I ) RE-P0063 

DO 2 J=2,JEND ReIp0064 

T *^ . RE-P0065 

Jlill. ^TOTNUM( II ) /PvJ RE-P0066 

"t*^"^ RE-P0067 

ri7 ifJ^ RE-P0068 

aj. FSTAT(D1,D2,ALPHA,F,IACC) RE-P0069 

RJKJhl = JNUII RE-P0070 

FFT7 = F 8((1.+ RINUI-I * PH0)/(1. + RJNUM * PHI(I2))) RE-P0071 

CA:,L PLEVEL(D1,D2,FP0W,P) RE-P0072 

IF(?.GE..9995D0) GO TO 12 ' RE-P0073 

PRINT 3, J.JNUM.FPOW.P RE-P0074 

FORMAT (IHO, ' NO. OF RANDOM LEVELS 15,' NO. OF PEOPLE IN A LEVRE-P0075 

*EL ,15,//' F VALUE TO CALCULATE POWER F8.3, ' POWER VALUE' ,F10RE-P0076 

*.5) 

GO TO 2 

PRINT 3, J.JNUM.FPOW.P 
GO TO 4 
CONTINUE 
CONTINUE 
CONTINUE 

RE-P00S4 
RE-P0085 



RE-P0077 
RE-P0078 
RE-P0079 
RE-P0080 
RE-POOCl 
RE-P0082 
RE-P0083 



ERIC 
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Appendix B-'-' 



^fSt^'^^^t^S^^I'^ °^ "'^ P'^'*^"" included here. The subroutines 
t t I^LEVEL are very long and were' excluded. FSTAT returns 

(F> gi^en alpha (ALPHA) degrees of freedom (Dl, D2) ; 
PLEVEL returns the probability of a given F statistic (P) given 
degrees of freedom (Dl. D2). Similar routim- should be available 
at most computer centers, 

2This program requires only one data card which contains: 

Columns Variable Format 

1~3 Alpha level f3.2 

0 F4.2 
8-11 F4.2 
12-15 RN, Sample size F4.0 
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C THIS 7/XOGRAM FINDS POWER VALUES FOR A TWO-WAY RANDOM-EFFECTS 

C ANALYfIS OF VARIANCE. 

C 

c 

PrUILE PR^qiSION D1,D2,ALPHA,F,FP0W,P 
R'iAD liAEPHA.PHI.PHO.RN 

1 FOEMT(F3.2,2F4.2,F4.0) 

VRiNT 2,alpi:a,phi,pho,rn 

2 FORMATCITHE VALUE OF ALPHA IS'F5.2//' THE VALUE OF PHI IS'F6.2//' 
COTHE VALUE OF PHI(SUBO), THE IiyPOIKEISI VALUE IS'F6. 2// 'OTHE KUMBE 
CR OF OBSERVATIONS IS'F6.0/////) 

N = RN 

KEND « RN/4.0 

lACC = 5 

DO 3 LLA = 1,3 

GO TO (20, 21, 22), LLA 

20 PRINT 6 

6 FORMAT ClTHE FOLLOIWING ARE POWER VALUES FOR THE INTERACTION'////) 
GO TO 16 

21 PRINT 7 

7 FORMAT ClTHE FOLLOWING ARE POVJER VALUES FOR FACTOR A'//' FACTOR A 
1 HAS I LEVELS'////) 

GO TO 16 

22 PRINT 8 

8 FORMAT ClTHE FOLLOWING ARE POWER VALUES FOR FACTOR B'// 
. 1' FACTOR B HAS J LEVELS'////) 

16 PRINT 1111 

1111 FORMAT (7X,'K, THE NO. 'T29, ' J,THE NO. Of',T49,'THE NO. OF' 
CT69,' NUMERATOR ',T86, 'DEN0MINAT0R'T104, 'F' ,T122, 'POWER') 

PRINT 1112 

1112 FORMAT (7X,'0F OBSERVATIONS' ,T29, 'LEVELS IN' ,T49, 'LEVELS IN' 
CT69, 'DEGREES OF' ,T86, 'DEGREES OF' ,T104, 'STATISTIC' ) 

PRINT 1113 

1113 FORMAT (7X,'IN A CELL' ,T29, 'FACTOR B' ,T49'FACTOR A' 
CT69, 'FREEDOM' ,T86, 'FREEDOM' ) 

DO 3 K = 2, KEND 
RK=K 

JEND - N/(K*2) 
DO. 3 J=2; JEND 
I '«.N/(K*J) 

GO TO (10, 11, 12), LLA 

10 Dl - (I-l * (J-l) 
D2 = i*j*(K-l) 

CALL FSTAT(D1,D2,ALPHA,F,IACC) 

FPOW « F * ((1. +RK * PH0)/(1. + RK*^''I)) 

GO TO 5 

11 Dl = I-l 

D2 - (I-l * (J-l) 

RJ = J 

CALL FSTAT(D1,D2,ALPHA,F,IAC(;) 
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FPOW = F * ((1. + RJ*RK*PH0)/(1. + RJARK*PHI)) 
GO TO 5 
12 Dl «=J-1 

D2 « (I-l) * 
RI = I 

CALL FSTAT(D1,D2,ALPHA,F,IACC) 
FPOW = F * ((1. + RI*RK*PH0)/(1. + RI*RK*PHI)) 
5 CALL PLEVEL(D1,D2,FP0W,P) 
PRINT 4,.K,J,I,D1,D2,F.P 

3 CONTUSE 
STOP 
END 



